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Big bang nucleosynthesis (BBN) and the cosmic microwave background (CMB) are two major 
pillars of cosmology. Standard BBN accurately predicts the primordial light element abundances 
( 4 He, D, 3 He and 7 Li), depending on one parameter, the baryon density. Light element observations 
are used as a baryometer. The CMB anisotropics also contain information about the content of the 
universe which allows an important consistency check on the Big Bang model. In addition CMB 
observations now have sufficient accuracy to not only determine the total baryon density, but also 
resolve its principal constituents, H and 4 He. We present a global analysis of all recent CMB data, 
with special emphasis on the concordance with BBN theory and light element observations. We find 
fl B h 2 = 0.02501^026 and y p = °- 250 ±o.oi4 (fraction of baryon mass as 4 He) using CMB data alone, 
in agreement with 4 He abundance observations. The determination of Y p allows us to constrain the 
relativistic degrees of freedom during BBN, measured through the effective number of light neutrino 
species, N„ >e ff = 3.02^0 79, in accord with the Standard Model of Particle physics. With this 
concordance established we show that the inclusion of standard, N Vif ,ff = 3, BBN theory priors 
significantly reduces the volume of parameter space. In this case, we find Qsh 2 = 0.0245lo'oo28 an d 
Y p — 0.2493^q'oo5q . We also find that the inclusion of deuterium abundance observations reduces 
the Y p and Jls/i 2 ranges by a factor of ~2. Further light element observations and CMB anisotropy 
experiments will refine this concordance and sharpen BBN and the CMB as tools for precision 
cosmology. 



I. INTRODUCTION 

Big-bang nucleosynthesis (BBN) and the cosmic mi- 
crowave background anisotropy (CMB) are two pillars of 
the hot big bang model. 

The theory of BBN has long stood as an emblem of 
the predictive power of the Big Bang model 0, 0, IE • 
BBN has long provided the most reliable measurement of 
the cosmological baryon density. However, the CMB is 
rapidly becoming the preferred method for determining 
the baryon density, with its rapidly increasing precision. 
With the CMB, light element abundance observations 
become a powerful probe of the universe. 

The observation and analysis of cosmic microwave 
background (CMB) anisotropies have attracted a great 
deal of attention in recent years due to their unique rele- 
vance for cosmological theory (see [j| for a recent review) . 
A flood of observational results have been published dur- 
ing the last two years H H H EE El El 13 These 
observations taken together measure CMB anisotropies 
over a large range of angular scales. The CMB is sensitive 
to the properties of the photon-baryon fluid and hence al- 
lows a precision determination of the baryon density at 
redshift z ~ 1000. 

It is therefore apparent that combining BBN and 
CMB provides an opportunity for meaningful consistency 
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checks on the standard cosmology and has the potential 
to be powerfully predictive probe of nuclear and par ticle 
astrophysics at low and high redshift EE EH EE HE EE 
19, 20]. Given the constraint on the baryon density from 
the CMB, BBN yields a tight prediction of primordial 
4 He abundance. 

We also explore the promise of combining the CMB 
data with measurements of the deuterium (D/H) abun- 
dance, showing that current measurements of (D/H) can 
combine with the CMB constraints to reduce the error 
bars on 4 He by another factor of 2. 

In section 2 we explain the data (CMB, BBN) and the 
method used to determine the likelihood surface in pa- 
rameter space - the Metropolis-Hastings Markov Chain 
Monte Carlo (MCMC) algorithm. In section 3 the pa- 
rameter confidence intervals that are extracted from that 
likelihood surface are discussed. In section 4 we discuss 
the implications of our results and how precision might 
be increased with further cosmological data. 



II. METHODS 

Our cosmological parameters are measured by use 
of the Markov Chain Metropolis Algorithm. We al- 
low the 4 He mass fraction to float as an indepen- 
dent variable, yielding the following parameter space: 
n m ,n Al n s ,h,n B h 2 ,Yp,T,r,n t where fi( m ,A,B) is the 
matter, cosmological constant, and baryon contents, re- 
spectively, Y p is the fraction of the baryons in 4 He by 
mass, h is the Hubble constant in units of lOOfcm/ sMpc, 
n( Sj t) are the power-law index of the primordial scalar 
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FIG. 1: The figure shows the data used in this paper (see legend) and 68% confidence intervals in the space of inferred theoretical 
power spectra (solid). 



and tensor perturbations respectively, r is the fraction 
of the observed CMB quadrupole that is tensor and r 
is the optical depth to the last-scattering surface (that 
is, re-scattering of CMB photons by reionization is al- 
lowed for). Note that in our selection of the parame- 
ter space simplifying assumptions have been made: adi- 
abatic, scale-free primordial perturbations, the universe 
contains only cold dark matter and a cosmological con- 
stant and the neutrino species are strictly those of the 
Standard Model. To determine the region of this parame- 
ter space allowed by experimental data, one must sample 
the space over a wide range of points. At a given point the 
relative likelihood of the parameter values yielding the 
observations must be determined, and the range of points 
sampled must adequately cover the space. Our primary 
likelihood calculation is a comparison of the simulated 
CMB spectra produced by CMBFAST against the 
WMAP CMB experiment 0, along with smaller-scale 
(bin i e ff > 350) data from the following experiments: 
Toco98 6], DASI 0, Maxima @L VS A @, ACBAR [l^, 
Boomerang02 HJ, and CBI [U. We include the 
published calibration uncertainties for each experiment 
and find the maximum likelihood value for these param- 
eters at each point in the cosmological parameter space. 

Figure 1 shows the data we used as well as the most 
likely 68% of our inferred theory power spectra. 

One could attempt to sample the parameter space 
on a uniform grid, but the high dimensionality coupled 
with the computational demands of CMBFAST makes 
this impossible in a reasonable amount of time. In- 
stead, we imp lemented the Metropolis-Hastings MCMC 
Algorithm [22I I23I |24| : starting from the current point 



in parameter space A*.; one proposes a test point 
drawn from a distribution described by a density func- 
tion called the proposal density: P{l?\ Aj). The choice 
of proposal density is somewhat arbitrary, but a poor 
choice will cause the parameter estimation procedure to 
be inefficient. The likelihood relative to A; is computed 
at the test point ~f?. If the likelihood at jf is greater, 
then point Xi+i = ~p\ Otherwise there is a probabil- 
ity that A*i + i = J? equal to the likelihood of ;? di- 
vided by the likelihood of AV If the proposal density 
is not symmetric in its arguments then this probability 
needs to be corrected by the factor P(X*i\ J?)/P(j?\ A*;) 
in order to enforce detailed balance. If A^+i is not set 
to the the point then it is set to the current point: 
Aj + i = A* \. After a sufficient number of iterations the 

resulting density of points | A^ | i — 1 . . . n| asymptot- 
ically approaches the likelihood function on parameter 
space. What constitutes a "sufficient" number of points 
is in general difficult to determine (and impossible to de- 
termine with absolute certainty) - though there are tests 
which are good indicators. We use a test for our Markov 
Chain which was suggested by [25| and also used by the 
WMAP team[26|. 

We tune our Markov Chain code for optimal efficiency 
by first finding the maximal likelihood point in parameter 
space using a global maximization method (within our 
prior space, see below) then using it as the starting point 
for a sample Markov Chain. The variance of the sample 
chain is used to compute a step size matrix that will 
be used by the main chains. An efficient Markov Chain 
should take steps that are not too large or too small - 
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either will result in an inefficient, slowly converging chain. 

We choose our proposal density to be a multivariate 
Gaussian, whose covariance matrix V s is proportional to 
the covariance of the sample chain V a . The proportion- 
ality constant a(D) is chosen so if the underlying distri- 
bution were Gaussian, 50% of the Markov chain points 
would be accepted. Its value depends on the number of 
dimensions of parameter space, D (for us D — 9). We 
found: 

a (D) ~ 0.54784 D - 0.36159 (1) 

We ran two sets of Markov Chains, A and B, each 
consisting of 20 independent chains. Each chain was 
started at a point chosen from the distribution of the 
sample chain. Set A had only very weak top-hat priors 
(n B h 2 G [0.014,0.030], Y p G [0.13,0.34], h G [0.45,0.95], 
fl m G [0.03,1.00], G [0.00,0.97], n s G [0.5,1.4], 
n t G [-3.0,3.0], r G [0.0,3.0], r G [0.0,1.0]), whereas set 
B additionally had two strong priors: a BBN consistency 
condition between fish 2 and K, |27| . and the constraint 
h = 0.74 + 0.11 - 0.094 gg. The BBN consistency 
condition is simply this: for a given baryon density one 
expects BBN to produce a certain abundance of 4 He, 
with some theoretical error (mostly driven by uncertain- 
ties in the nuclear cross-sections). In set A we treat Q,sh 2 
and Y p as two independent variables. In set B we enforce 
theoretical self-consistency between those variables. This 
non-CMB constraint is incorporated in a Bayesian way. 
The additional information is included as an additional 
prior and the density to be sampled from (the posterior 
density) becomes 

P tot (?) = Pcmb (?) Pbbn (tt B h 2 , Yp) P Hubble (ft) 
= P (j£\CMB^ P (Y p \n B h 2 ) P {h\Hubble data) 

(2) 

where P (x\y) is the conditional density for getting x, 
given y. Note that P (Y p \VlBh 2 ) is purely a theoreti- 
cal prior enforcing the BBN relation, and contains no 
abundance measurements. The Markov Chain then au- 
tomatically explores the new, more constrained region of 
parameter space. As one might expect, sets A and B dif- 
fer significantly in their parameter space coverage, and 
thus their proposal densities and chain starting points 
were determined separately. 

Though a Markov Chain approach saves significant 
computational time, it is difficult to guarantee after some 
number of points that the chain has converged sufficiently 
to the true, underlying distribution. Indeed, a chain can- 
not tell one anything about a region it has not yet visited. 
We use a convergence test suggested by Gelman & Ru- 
bin , which was also employed by WM AP |2(| . 

We have generalized this criterion to multiple dimen- 
sions, keeping in mind that any convergence test must 
be covariant (if a transformation of parameter space can 
change the determination of "convergence" , then the test 
is a bad one) . Each chain out of the 20 has its own mean 
and variance . If each chain reflected the underlying 



distribution, then the variance of the means of the chains 
should be much less than the variance of the underlying 
distribution. We thus compute the variance of the chain 
means, multiply with the inverse of the average chain 
variance, and take the trace: 




3=1 



N H] 
j = l J i=l 

H = Tr [UW" 1 ] jD 

(3) 

where N is the number of chains (20), rij is the number 

of points in chain j, X is the total mean, X j is the mean 
of chain j. We require that [i < 0.1. Set B easily satisfies 
this criteria with 30, 000 points, whereas set A required 
about 60,000 points. The average chain variance, W, 
is used because this underestimates the variance of the 
distribution until convergence is attained. 

As a self-consistency check, one can take the point dis- 
tribution of set A and combine it with the BBN fish 2 — Y p 
and Hubble Key project priors. The resulting distribu- 
tion should be the same as set B. The extent to which 
these distributions differ is a measure of non-convergence 
of the sets. We determined that the 68% confidence re- 
gions of these distributions more than 95% overlap in 
the (yish 2 , Y p ) plane. A Markov Chain can be combined 
with a prior after the generation of the chain by assign- 
ing a weight to each point. The likelihood of a region 
in parameter space is then the weighted density of the 
points in that region. Because a Markov Chain main- 
tains the full D-dimensional likelihood distribution in 
the parameter space, after it is generated the chain may 
be convolved with any arbitrary other likelihood func- 
tion in that parameter space. Thus one can generate a 
Markov Chain distribution for WMAP alone, and chose 
any subset of the other cosmological datasets to convolve 
it with - for very little additional CPU cost. This is the 
basis of the Cosmic Concordance Project (CCP) web- 
site http://galadriel.astro.uiuc.edu/ccp/ where 
the parameter constraints from chains A and B can be 
explored and combined with other cosmological datasets 
and priors. Further details about the CCP and about 
our parameter estimation methodology will be given in 



III. RESULTS 

As our Markov Chain sets A and B approached 60, 000 
and 30, 000 points respectively, convergence test eq. 
gave fi ^ 0.05 — 0.06 01 an d we declared our chains 
sufficiently converged to provide reliable statistics. Fig- 
ure 2a shows the 68% confidence region of set A in the 
{fish 2 ,Y p ) plane. Thus from CMB data alone we find 
n B h 2 = 0.025^S and Y p = 0.250i°;°i° 
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FIG. 2: The top figure shows the 68% confidence region in the Qph 2 , Y p plane from the CMB data alone (WMAP and high-£ 
data from other recent CMB experiments) bounded by the blue line. The black line bounds the 68% confidence region from 
BBN theory alone, using the best fit Q,ph 2 from the CMB. Note that the BBN theory band is in good agreement with the 
CMB data. The bottom figure shows the result of combining CMB and BBN data using an expanded Yp axis. The solid red 
region is the 68% confidence region for the set B Markov chains which have as priors the BBN constraint and the Hubble Key 
project constraint on h. Note that the BBN constraint greatly reduces the allow range of Yp as a function of fis/i 2 . 



It is worth noting also that Markov chain set A yields a 
tight constraint on the neutrino number during BBN: Al- 
lowing the number of BBN neutrinos to float, one can find 
a 68% confidence interval for the number of BBN neu- 
trinos needed to yield the determined Y p from the deter- 
mined fish 2 . Given a BBN code that computes the prob- 
ability density P(Y p \Vl B h 2 , N u ^ e ff) (where the stochastic- 
ity is due to the measurement uncertainty in the relevant 
nuclear cross-sections) we can compute P(N u>e ff) based 
on CMB data and BBN theory as 

P(A^ e// |BBN,CMB) = 

J dY p d(n B h 2 )p(N„, eff \Y pi n B h 2 )p(Y pi n B h 2 ) = 
J dY p d(n B h 2 )p(Y p \n B h 2 ,N^ eff )p{Y p ,n B h 2 ). 

The first term under each integral enforces the 
BBN relation and the second term the CMB pos- 



terior. The second equality holds true if we as- 
sume flat priors P(Y p \il B h 2 , BBN) — const, and 
P{N Utf ,ff\Vt B h 2 ,BBN) = const for the range of param- 
eter space of interest 0- We nnd N »,eff = 3-0 2 tol?9- 
N v>e ff is consistent with the standard model value of 3 
and previous studies E E E E E El ■ 

However, fl B h 2 and Yp are jointly constrained by BBN 
theory, and thus are not really independent variables. 
Adopting the standard BBN model (AT„, e// (BBN)=3) 
of [27| yields a consistency relation between Q B h 2 and 
Yp. The 68% confidence region of this consistency re- 
lation appears in Figure 2a as a narrow band (narrow 
enough that the upper and lower bounding curves ap- 
pear to merge). Enforcing this condition greatly in- 
creases the precision of parameter estimation, as evident 
in table 1, with dramatic affect on Yp measurement: 
Y P = 0.2493±gX. This is simply a result of the ac- 
curacy with which 4 He is determined by BBN (~ 0.1%). 
In Figure 2b we have zoomed in on this CMB-BBN con- 
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FIG. 3: The primordial abundance of deuterium can be used to further constrain the baryon density (blue shaded band is 
68%), which in turn increases the precision of the Yp determination (solid red region is 68%). Two different D abundances are 
used: 2a) data from 2 multiple-line absorption systems, and 2b) data from 5 systems (the 2 multiple-line plus 3 others) (see 
text). 



cordance region. Also shown is the 68% confidence re- 
gion of set B shaded as solid red. As one would expect, 
set B agrees with the product of the CMB (set A) and 
BBN (consistency band) likelihoods. It is important to 
note that in Figure 2a the agreement between the CMB 
and BBN allowed regions need not have happened. In- 
stead, the BBN consistency band might not have passed 
through the high CMB likelihood region, which would 
have forced one to consider a BBN scenario other than 
the standard model one with 3 neutrinos. The CMB- 
BBN agreement reaffirms the standard BBN scenario. 

Model A and Model B compare quite well to 4 He ob- 
servations. Olive, Skillman and Steigman (1997) [3(j and 
Fields and Olive (1998) [3l| find Y p = 0.238±0.002, while 
Izotov and Thuan (1998) [H find Y p = 0.244 ± 0.002. 
The errors cited are statistical only. Comparing these 
numbers, not only are they discrepant from each other, 
but they lie below the mean value determined in this 
evaluation. However, Olive and Skillman (2001) crit- 
ically evaluate the methods used in determining Y p and 
find a lower bound to a systematic error, <j sys > 0.005. 



This systematic error is added in quadrature with the 
statistical error to determine the total error, increasing 
the errors to 0.0054. The Y p observations are brought 
into marginal accord with each other and the CMB; both 
lie systematically lower than the CMB determined value. 
As discussed earlier, the systematic error used is only a 
lower bound, and as such the true errors are most likely 
larger than those quoted. 

In Figure 2 the allowed fi^/i 2 is very large. Any other 
data that can reduce the allowed VLb h 2 range will have 
the additional benefit of refining the precision of the 
Yp measurement. As an example, we consider D abun- 
dance in Figure 3. The value of Dj H is still a somewhat 
open question due to low number statistics, and thus we 
demonstrate the effects of two different D abundances. 
In Figure 3a we use the average of the 2 multiple absorp- 
tion line systems of H3 : D/H = (2.49 ± 0.18) xl0~ 5 . 
The blue shaded band is the 68% confidence baryon den- 
sity range allowed b y th is D/H value as determined with 
the BBN theory of 27]. The solid red region is the 68% 
confidence region resulted from the convolution of the D 
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TABLE I: Cosmological Parameter Estimates 
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Model A 
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0.0250_ 0026 


0.0245_ 0029 


V 


n ot;n+o.oio 
0.250_ 0014 


n o/inq+0.0007 
u.iiao_ 001(] 


n 


n R«/i+ ' 057 


U. Ma ± U.UOy 


ft M 


n 241 +0.062 
U.Z41_ 064 


n 21 q+ 041 

u.ziy_ 05g 


V A 


n 7Q2+ 063 


7Q8+ 060 

u - ' y °-0.045 


n s 


1.047lHf 5 


i n24 +0 041 


r 


0.168lg-°g 


0.117in? 4 7 


n t 


0.0891°;$ 


0.2071°;^ 


T 
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baryon range and the concordance region of Figure 2b. 
Thus we find Y p = 0.2491 4 ; ; 0004 , fl B h 2 = 0.0237±°' g * . 
Alternatively, one can use a conservative combination of 
5 D/H measurements, including the two multiple ab- 
sorption line systems employed above [34 , l35t 139 . 137L l38| : 
D/H = (2.78 ±0.29) x 1CT 5 (the overall error increase 
is because the other three systems are not consistent 
with each other or the multiple absorption line systems, 
a hint of an underlying systematic error for the single 
absorption line systems). For this D abundance we find 



Y p - U.2l?vS iU - 



0004 and n B h 2 = 0.02301°;°°^. 



IV. CONCLUSIONS 

This work has been based on two general ideas; (1) 
BBN and the CMB independently probe two different 
epochs, providing valuable consistency checks for the un- 
derpinnings of the Standard Cosmological Model; (2) 
having established that the cosmological model agrees 
remarkably well with these very different observational 
probes we use these data to to make a precision mea- 
surement of the 4 He abundance. We have presented an 
analysis of all recent CMB data, in which we have de- 
termined the cosmic baryon density and the primordial 
helium abundance. We found £l B h 2 = 0.02501^^6 and 

Y p = 0.250lHi° at 68% from CMB data alone - This is 
consistent with 3 standard model neutrinos during BBN. 

We have shown that this is fully consistent with 
the predicted 4 He abundance from BBN theory, and 



marginally consistent with 4 He observations. The likely 
source of this slight discrepancy is an underestimate of 
the dominant, systematic uncertainties in the 4 He obser- 
vations, which now seems affirmed with the CMB deter- 
mination of Y p . The agreement between the CMB-only 
set A confidence region of fl B h 2 , Y p and the consistency 
band based on BBN theory shown in Figure 2a reaffirms 
the standard BBN model. Thus using BBN theory, we 
can effectively remove Y p as a freely floating variable, en- 
forcing the n B h 2 -Y p BBN relation in CMB data analysis. 
Given this, we found the incorporation of BBN theory 
into parameter extraction from CMB data results in a 
precision measurement of the 4 He abundance. We find 
fl B h 2 = 0.02451°;°°^ and Yv = 0.24 93+ 0007 



-o.ooio- 



Using this CMB-BBN determined Y p , one can study 
and constrain stellar evolution [ImJ^cJ. One can similarly 
study the nucleosynthetic history of all the light element 
abundances as discussed in ^(| and references therein. 

We show the promise of incorporating deuterium 
abundance observations, yielding Y p — 0.249llg qq 04 , 
fl B h 2 = 0.0237l° 1 00 1 - or Y p = 0.2488 4 ; ; 0004 , Q B h 2 = 
0.0230tQ; 00 °2 depending on which systems are used to 
measure the deuterium abundance. 

The addition of the D abundance observations are only 
one example of many possible cosmological datasets that 
might be incorporated into parameter extraction to in- 
crease precision. An experiment may constrain a param- 
eter directly, or may reduce degeneracy in a related pa- 
rameter. For example, using large-scale structure infor- 
mation to reduce the residual n s -fl B h 2 degeneracy in the 
current CMB data would also increase the precision of 
the Yp determination. Also, further light element obser- 
vations and CMB anisotropy experiments will refine this 
concordance and sharpen BBN and the CMB as tools 
for precision cosmology. Due to the effect of the 4 He 
abundance on the damping tail, this may improve the 
constraint on a possible running of scalar spectral index. 
These are the topics of on-going work and can be further 
explored at the Cosmic Concordance Project web-site: 
http : / / galadr iel . astro . uiuc . edu/ccp/| |29l | . 
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Science Alliance under grant AST020003N and utilized 
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